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There is great interest in designing photonic devices capable of disorder-resistant transport and 
information processing. In this work we propose to exploit 3D integrated photonic circuits in 
order to realize 2D discrete-time quantum walks in a background synthetic gauge field. The gauge 
fields are generated by introducing the appropriate phase shifts between waveguides. Polarization- 
independent phase shifts lead to an Abelian or magnetic field, a case we describe in detail. We 
find that, in the disordered case, the magnetic field enhances transport due to the presence of 
topologically protected chiral edge states which do not localize. Polarization-dependent phase shifts 
lead to effective non-Abelian gauge fields, which could be adopted to realize Rashba-like quantum 
walks with spin-orbit coupling. Our work introduces a flexible platform for the experimental study of 
multi-particle quantum walks in the presence of synthetic gauge fields, which paves the way towards 
topologically robust transport of many-body states of photons. 


Introduction. A long-standing aim in condensed mat¬ 
ter physics is to understand the behavior electrons in 
two dimensional systems in the presence of a magnetic 
field [1]. The reasons for this are both of fundamen¬ 
tal and of applied nature. When the system is well de¬ 
scribed by weakly interacting quasielectrons, it is known 
that topologically-protected edge states akin to those of 
topological insulators [2] are present. Strongly interact¬ 
ing electrons in a magnetic field arrange themselves in 
non-standard states of matter [3] which cannot be de¬ 
scribed by a local order parameter and present topologi¬ 
cal order [4] . The excitations of this state of matter may 
present non-Abelian statistics, which could be used for 
topologically-protected quantum computation [5]. 

The promise of ground-breaking applications together 
with the richness of the underlying physics of two- 
dimensional (2 D) quantum particles in a magnetic field 
has made these systems a favorite subject of quantum 
simulator proposals [6]. In these quantum simulators— 
physical systems unnaturally made to behave according 
to a specific model —the magnetic field is artificial, i.e. 
synthetic. Instead of using charged particles in an ac¬ 
tual magnetic field, in a quantum simulator one typi¬ 
cally uses neutral particles upon which the effects of a 
fictitious magnetic field are imposed. For neutral cold 
atom approaches, methods used to generate a synthetic 
magnetic field include rapid rotation [7, 8] , Raman-laser- 
induced Berry phases [9], laser-stimulated tunneling in 
optical lattices [10-15], or lattice shaking [16]. 

An alternative approach to quantum simulation is to 
directly implement the time-evolution of the system, 
as opposed to engineering the underlying Hamiltonian. 
Quantum walks (QW) [17] are a prominent example of 
this idea and have been realized in a variety of plat¬ 
forms, including neutral trapped atoms [18], trapped ions 
[19, 20], and nuclear magnetic resonance in continuous 
[21] and discrete-time [22]. A promising platform are 


photonic quantum simulators [23] , which have been used 
to simulate QWs in the bulk [24, 25] and in waveguide lat¬ 
tices [26], as well as photon time-bin encoded QWs [27]. 
Furthermore, two-particle QWs [28] have been realised 
in integrated photonic circuits using quasi-planar geome¬ 
tries [29-31], non-planar circuits in a “criss-cross” con¬ 
figuration [32, 33], and Anderson localization has been 
reported in the disordered case [34] . 

Discrete-time QWs (DTQWs) in ID may be imple¬ 
mented with a planar integrated photonic circuit (IPC) 
forming an array of beam-splitters [31]. Each beam¬ 
splitter performs the coin and step operator at the same 
time, shifting the photon left and right in quantum su¬ 
perposition. Successive beam-splitters create further su¬ 
perpositions, leading to the genuinely quantum interfer¬ 
ence phenomena which are characteristic of QWs. In this 
implementation, time is encoded in the direction of prop¬ 
agation of the photon in the IPC. 

A promising development in IPC technology is the ca¬ 
pability to print the waveguides in a truly 3D configu¬ 
ration. In particular, it is possible to implement quan¬ 
tum walks on a 2D lattice using a 3D network of beam¬ 
splitters. In such a network, each waveguide corresponds 
to one lattice site, and there are vertical and horizon¬ 
tal beam-splitters, which shift the photon wave function 
in an up-down and left-right superposition, respectively 
(see la). Similarly to the implementation of ID quan¬ 
tum walks, the time is encoded in the spatial direction of 
propagation of the photon. 

In this work we propose 3D integrated circuits to re¬ 
alize 2D QWs in a synthetic gauge field. This is ac¬ 
complished by introducing controlled phase-shifts be¬ 
tween waveguides at the beam-splitters. The phase-shifts 
are chosen in such a way that the photons gain global 
phases when going around a closed loop, leading to the 
Aharonov-Bohm effect [35] (see Fig. (lb)). Polarization- 
independent phase shifts lead to an Abelian or magnetic 
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FIG. 1: (a): Diagram of the proposed 3 D photonic circuit 
that realizes the DTQW on a 2D lattice. The z axis rep¬ 
resents the direction of time, while the xy plane represents 
the two spatial dimensions where the QW takes place. Each 
full time-step of the quantum walk is implemented by a se¬ 
quence of four sub-steps, depicted here and in (c) in four 
different colors. Each substep implements a set of mutually 
commuting gates between waveguides which can be applied 
simultaneously. The couplings in the y direction introduce 
a ^-dependent phase-shift between waveguides which mimics 
the effect of a magnetic field, (b): The accumulated phase 
acquired by the quantum walker going around two examples 
of closed trajectories with opposite chirality. This phase de¬ 
pends only on the chirality and the number of elementary 
cells inside the loop, as in the Aharanov-Bohm effect expe¬ 
rienced by a charged particle in a constant magnetic field, 
(c): The four different kinds of links involved in the quantum 
walk. Each group of links is depicted in a different color. The 
DTQW is realized by creating superpositions between sites in 
each group in a sequential way. 


field, while polarization-dependent phase shifts lead to a 
non-Abelian gauge field [36] . Our scheme may be readily 
generalized to QWs of two or more photons [28], allowing 
for the implementation of QW exhibiting topological fea¬ 
tures [37, 38] in the multi-photon case in an IPC. Further¬ 
more, the spatial dependence of the effective gauge field 
is highly tuneable, thus allowing for synthetic gauge fields 
in exotic configurations, such as magnetic monopoles, 
with no added difficulties. There is great interest in en¬ 
gineering photonic technologies with topologically pro¬ 
tected properties [39] . Although several examples of pho¬ 
tonic systems with topologically protected edge states 
have been proposed [40-42] and realized [43, 44] with 
laser light, such as the quantum Hall effect and the Flo- 
quet topological insulator [45], our proposal is to realize 
quantum walks in effective gauge fields in the few-walker 
regime, using single-photons. 

2D QW in a synthetic magnetic field with an IPC. 
The evolution of a charged bosonic particle in a 2D lat¬ 


tice with a perpendicular magnetic field is described by 
the Hamiltonian 

H = J ^ a }n,n+l a m,n + a m+l,n a rn,n + H.C.) . ( 1 ) 

m,n 

The operators n and a m?n create and destroy one par¬ 
ticle at site (m,n) of the lattice, respectively, and obey 
bosonic commutation relations. The constant J is an 
arbitrary energy scale and is the magnetic flux per 
plaquette. The key feature of this Hamiltonian is that 
hopping in one of the directions of the lattice entails 
the acquisition of a position-dependent phase. The spe¬ 
cific spatial profile of these phases is such that the global 
phase acquired by a particle going around a closed path 
on the lattice is position- independent and equal to e l( ^ N , 
where N is the number of elementary cells inside the 
path (see Fig. (lb)). The particular choice of phases 
is arbitrary—in Eq. (1) we chose the so-called Landau 
gauge for convenience—as long as the accumulated phase 
along closed paths leads to the correct global phase. The 
idea of introducing position-dependent phases has previ¬ 
ously been used to realize chiral QWs on graphs [46, 47], 
as well as a QW in an effective electric field [48, 49]. 

Here, we use an approach involving coinless discrete¬ 
time quantum walks on a 2D lattice where each step im¬ 
plements a position-dependent phase, in analogy with the 
dynamics given by Hamiltonian from Eq. (1). We ex¬ 
plain how to implement this quantum walk in a 3D IPC 
and present numerical evidence that its dynamics shows 
similar features to the one described by (1), namely the 
presence of topologically protected edge states. Although 
in other photonic implementations of discrete-time quan¬ 
tum walks the polarization of the photon is used as the 
coin [25, 27], here we assume the IPC to be polarization 
independent so that we can use entanglement in the po¬ 
larization to simulate bosonic and fermionic statistics [28], 
as previously done in ID quantum walks in Refs. [34], and 
so we have a coinless quantum walk. For a general def¬ 
inition and discussion of spreading properties of coinless 
quantum walks on lattices see Ref. [50]. 

To implement the quantum walk on the 2D lattice us¬ 
ing a 3D IPC, a lattice site in a position (x,y) will cor¬ 
respond to a waveguide engraved in the IPC, also la¬ 
belled by the position (x,y), which is extended in the 
z-direction, corresponding to the time dimension of the 
DTQW. The hops of the quantum walk correspond to 
sequences of beam splitters, which can be implemented 
by bringing adjacent waveguides close together in such 
a way that they are evanescently coupled. Each lattice 
site has at most 4 neighbors but, since we can only cou¬ 
ple two waveguides at a time, it is not possible for the 
photon in a certain waveguide to hop to all its nearest 
neighbors in one step. This way, we divide the DTQW 
into 4 substeps, as depicted in Fig. lc, where the links 
in green correspond to hopping terms that will be imple¬ 
mented by a unitary matrix C/i, the ones in yellow by C/ 2 , 
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the ones in red by Us and finally the ones in blue by C/ 4 . 
This way, one full step of the quantum walk will be given 
by 

b^step U A U 3 U 2 U 1 . ( 2 ) 

This sequence of unitaries can be applied many times 
along the z-direction of the IPC to implement subsequent 
steps of the quantum walk. A schematic representation 
of part of the 3D photonic circuit implementing the red, 
green and yellow links is shown in Fig. la. 

In order to mimic the effect of a magnetic field, we need 
to construct U step in such a way that if the walker makes 
a closed path around N elementary cells, it acquires a 
position-independent phase e lc ^ N . To show the form of 
the unitaries U{ which satisfy this requirement, we define 
the states of the Hilbert space \x)\y), with x G {1, M} 
and y G {1, M}, and say that the photon is in state \x)\y) 
if it is in the waveguide labelled by the coordinates (x,y). 
In this basis, we define the hopping operators in the x 
and y directions as 

Vx=-^(\x){x\ + \x+l){x + l\) (3) 

+ ^(|z)(z + 1 | + |z + 1 )(:z|) 

V v(<t>)=-j=(,\ y ){ y \ + \y + i)(y + M) ( 4 ) 

+ -L(e- i <t>\y)(y+l\ + e i *\y+l)(y\), 

which corresponds to an unbiased beam-splitter matrix 
and to a phase-shifted beam-splitter, respectively. The 
operators J7i, C/ 2 , C /3 and C /4 are then defined as 

M/2-1 M/2-1 

C /1 = V2x+1 C /2 = V 2x 

x=0 x=l 

M M/2-1 

u 3=Y, E k>(*l®*W^) (5) 

x=l y= 0 
M M/ 2-1 

^=£ Tl \x){x\® V 2y{x4>). 

x=l y=1 

Us and C /4 cause the hopping of the photon in the y- 
direction and apply a phase which is proportional to 
the coordinate x. Previous experiments have shown full 
phase-shift controllability between two waveguides, by 
deforming one of the waveguides and thus creating a dif¬ 
ference in the optical path length [34]. Hence, the exper¬ 
imental implementation of C/ Step , although challenging, is 
within reach of current technology. 

Signatures of the magnetic field for single photons. To 
confirm qualitatively the correspondence between the 
proposed IPC and the time evolution generated by 
Eq. (1), we have computed the spectrum of the effec¬ 
tive Hamiltonian H e ff = ilogC/ s t ep , where U s t ep is the 
unitary operator implemented by the proposed optical 
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FIG. 2: In an ordered lattice, a magnetic field hinders the 
spreading of the wave function and quantum transport; in a 
disordered lattice, a magnetic field enhances transport due to 
the presence of non-localized edge states. The variance of the 
single-photon wave function vs. number of time-steps, for dif¬ 
ferent values of the magnetic flux <f is displayed in (a). The 
inital state is localized at the center of the lattice. For non¬ 
zero <f the QW spreads at a slower rate. In (b), transport 
efficiency 77 as function of time in the presence of Gaussian 
disorder with strength 5 = 0.1 (defined in the Appendix). 
The quantity 77 measures the accumulated norm at the oppo¬ 
site corner, or the efficiency of quantum transport across the 
lattice. In the presence of disorder, the applied magnetic field 
qualitatively enhances transport. 

circuit. In the Appendix we plot the spectrum of H e ff 
as a function of <f . We obtain a figure very similar to 
Hofstadter’s butterfly [51]. 

Furthermore, we investigate the effect of the synthetic 
magnetic field on the spreading and transport properties 
of the QW at the single-photon level, with and without 
disorder. Controlled disorder may be implemented via 
small, random differences in waveguide lengths at the 
evanescent couplings [34], which lead to fluctuations in 
each waveguide’s optical path (see Appendix). These 
fluctuations are static, in the sense that they are not 
time-dependent (^-dependent in Fig. (la)). 

To determine how quickly the QW spreads without 
disorder, in Fig. (2a) we plot the variance of the single¬ 
particle probability distribution a 2 as a function of the 
number of steps, for different values of 0 , the magnetic 
flux per plaquette. The initial wave function is localized 
at the center of the lattice. Although we plot here the 
result for three values of 0 , we have observed that the 
variance is always smaller for 7 ^ 0 than for </> = 0 . 
Hence, without disorder the magnetic field is detrimen¬ 
tal to the expansion of an initially localized photon wave 
function. We also study the QW evolution by comput¬ 
ing the transport efficiency between two far-apart waveg¬ 
uides in the presence of disorder. We choose one corner 
of the lattice as an initial site/waveguide, and the site at 
the opposite corner as the target. We introduce absorp¬ 
tion at the target waveguide, corresponding to position 
(M, M), at each step of the QW by replacing the opera¬ 
tor U step by U step (X — \M)\M)(M\(M\). Our measure of 
transport efficiency, 77 , is the accumulated probability of 
finding the photon at the target, 77 = ^ 2 t |(target |t/(£))| 2 , 
an approach similar to that used in [52]. Let us stress 
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that rj could be measured in an experiment, by coupling 
the target waveguide to a long chain of waveguides as 
proposed in [33], or to a detector at every time step. 

With disorder, low transport efficiency is expected, 
due to Anderson localization. Interestingly, we find that 
while in the ordered case a magnetic field slows down 
the expansion of the QW, in the case of disorder it does 
the opposite, thus enhancing quantum transport (see 
Fig. (2b)). This is attributable to the presence of chi¬ 
ral edge states (see Fig. (3) in Appendix ) which are 
topologically protected against localization. 

2-photon QW with a magnetic field. The single-particle 
probabilities obtained by using one photon as the input 
state of the IPC can be reproduced by using a classi¬ 
cal laser light source. However, if two or more indis¬ 
tinguishable photons are used as input, the probability 
distribution measured at the output of the circuit has 
no classical analogue and for many photons is, in gen¬ 
eral, hard to calculate [53, 54]. Also, by choosing appro¬ 
priate entangled states of two-photons, the statistics of 
bosons and fermions can be mimicked [28] and bunch¬ 
ing and anti-bunching phenomena have been observed 
in 1-D QW’s [34]. Here, we compute some observables 
for the DTQW’s of two entangled photons in a synthetic 
magnetic field. The average distance between photons is 
plotted in Fig. 3 for two entangled photons starting at 
the corner of the lattice. It is clear the effect of the par¬ 
ticle statistics in this quantity since bosons remain closer 
than fermions. Also, the presence of the magnetic field 
increases the average distance between particles. The 
presence of two-particle edge states can be seen from the 
probability that both photons are at the edge of the lat¬ 
tice shown in the Appendix. 

Non-Abelian 2D QW. Our proposed scheme to real¬ 
ize a magnetic QW with IPC may be generalized to a 
non-Abelian magnetic QW, provided the relative phases 
between adjacent waveguides are made polarization- 
dependent in a controlled way. When polarization is 
taken into account, a general term coupling two adja¬ 
cent lattice sites i and j can be written in the form 
a\^Jdj T + H.c. , where £, r run over photon po¬ 
larizations, and now a\^ creates a photon in site i with 
polarization £. 

Thus, to realize a QW in a non-Abelian synthetic gauge 
field, the beam-splitter matrices must be polarization- 
dependent, which are now described by 4 x 4 matrices 
instead of 2 x 2. In general, the beam-splitter matri¬ 
ces corresponding to different links of the lattice will not 
commute with each other, and will lead to non-trivial 
non-Abelian fluxes when the photons go around a closed 
loop (see Fig. (lb)). This is tantamount to a modified 
Aharonov-Bohm effect, where the photon wave function 
is multiplied by the Wilson loop [55] instead of a phase. 

Remarkably, interesting non-Abelian QWs may be im¬ 
plemented using relatively simple IPCs. In particular, a 
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FIG. 3: Average distance between photons for the two photon 
quantum walk on the IPC, after 20 steps, in a lattice of size 
30x30, for different values of the magnetic flux 0. The ini¬ 
tial state is localized at positions (1,1) and (1,2) of the lattice 
and the photons’ polarization states are entangled in a sym¬ 
metric (antisymmetric) way so that the exchange statistics of 
the wave-function is bosonic (fermionic). We observe, as ex¬ 
pected, that for bosonic statistics the photons tend to remain 
closer than in the fermionic case. Also, the presence of the 
magnetic field increases the distance between the photons. 

QW with Rashba spin-orbit coupling [56] may be realized 
with the choice U x = exp (iacr y ) and U y = exp(— iacr x ), 
where a x and a y are Pauli matrices. In the 3D IPC archi¬ 
tecture, this means adjacent waveguides in the x direc¬ 
tion are coupled with U x and those in the y direction with 
U y . Note that this choice does not require position depen¬ 
dent delays between waveguides as the Abelian magnetic 
field case does. In this scenario, since the circuit is now 
polarization-dependent, it would not be possible to simu¬ 
late different particle-statistics by entangling photons in 
polarization. 

Conclusion. We have introduced a scheme that allows 
implementing quantum walks in synthetic gauge fields 
using integrated photonic circuits. This scheme requires 
a strong experimental and technological effort: we need 
the capability to engineer 3D structures with a signifi¬ 
cant number of steps. In the last year several improve¬ 
ments have been achieved: 8 mode fast Fourier transform 
with 3D structure [57], reconfigurable phase [58] and op¬ 
eration at telecom wavelength which ensures lower losses 
and hence the possibility to realize longer chips [58] . Our 
proposal is well suited for the study of topological insu¬ 
lators at the single and few photon levels and it is highly 
flexible, allowing for the simulation of both Abelian and 
non-Abelian gauge fields. We have studied the single¬ 
photon quantum walk in a constant Abelian or mag¬ 
netic field and computed experimentally-accessible ob¬ 
servables, demonstrating topological properties, namely 
the presence of edge states enhancing transport across 
disordered lattices. We have also computed observables 
for two-particle quantum walks that demonstrate the role 
of entanglement and magnetic field in the behaviour of 
the walk. Overall, we have shown that the development 
of 3D integrated photonics can lead to the experimental 
study of interesting 2D quantum physics with topological 
features in the few-body regime. 
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APPENDIX 

Spectrum of the unitary implemented by the Integrated Photonic Circuit 

We have constructed a unitary matrix which implements one time step of the discrete-time quantum walk (DTQW) 
with a synthetic magnetic field, defined by Eq. (2) of the main text. This matrix, denoted as U step , can be decomposed 
in a product of beam-splitter and phase shifter matrices, which can be implemented in an Integrated Photonic Circuit 
(IPC). In 4 we plot the spectrum of the effective Hamiltonian H e ff = ilogU step as a function of <f>. As expected, we 
obtain a figure very similar to Hofstadter’s butterfly [51] —a complex, self-similar structure which arises in the case of 
electrons propagating on a 2D lattice in a strong magnetic field. For relatively large lattices (30 x 30), the spectrum 
of H e ff presents a fractal nature and a structure of gaps which is very reminiscent of Hofstadter’s butterfly. 
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FIG. 4: Spectrum of the effective Hamiltonian H e ff = zlogUstep which generates the Abelian QW proposed here, as a function 

of the magnetic flux per plaquette, f>. 



Time-evolution of single-particle quantum walks in a magnetic field 

We study the effect of a synthetic Abelian gauge field on the transport properties of a discrete-time 2D QW in the 
presence of disorder. All simulations have been done with the beam-splitter matrices given in Eq. (4) and Eq. (5) 
of the main paper. The initial wave function is localized at the lower left corner of the lattice. In Fig. (5) we plot 
the evolution of the QW without magnetic field and without disorder. There is efficient transport from one corner of 
the lattice to the other. In Fig. (6), localization close to the initial position is observed for zero magnetic field and 
disorder; transport from one corner of the lattice to the other is highly inefficient in this case. In Fig. (7) we plot the 
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evolution of the QW for non-zero magnetic field, </> = J, and in the presence of disorder. Transport from one corner 
of the lattice to the opposite one is clearly accomplished by edge states, which do not penetrate significantly into the 
bulk of the lattice. 



FIG. 5: Simulation of the IPC implementation of a DTQW in a 2D lattice in zero magnetic field without disorder, starting at 
position (X,Y) = (1,1). Probability distribution of a discrete-time QW on a 60x60 ordered lattice after 20, 40 and 60 time steps. 
The quantum walk propagates quickly from one corner of the lattice to the opposite corner. The white color corresponds to 
probabilities above 0.01, which can go up to 0.11 in (a), 0.056 in (b) and 0.016 in (c). Total probability for each plot: (a) 1 , 

(b) 1, and (c) 0.77. 




FIG. 6: Simulation of the IPC implementation of a DTQW in a 2D lattice in zero magnetic field with disorder, starting at 
position (X,Y) = (1,1). Probability distribution of a discrete-time QW on a 60x60 disordered lattice after 20, 40 and 60 steps. 
The strength of disorder is a = 0.2 and the magnetic field is set to zero. The presence of disorder hinders the propagation on 
the lattice and the photon remains localized close to its starting position, resulting in a low transport efficiency to the opposite 
corner. The probability amplitude is averaged over 20 disorder realizations. Note that the scales of the plots are different, for 
better visualization. Total probability for each plot: (a) 1, (b) 1, and (c) 0.997. 

At each step of the QW, the probability at the target waveguide is subtracted from the wave function, as we 
introduce absorption by replacing the operator U with U exp(—al arget a target ). Thus, the deviation from unity of the 
total probability is equal to the transport efficiency, 77 = | (target \^{t))\ 2 . 

The introduction of static disorder in the DTQW in a 2D lattice is done by multiplying each beam-splitter matrix 
V x or V y (0), defined in Eqs. (3) and (4) of the main text, by a matrix with random phases in the diagonal, 


V 

y o 


dis. 

x,y 




v* 


x,y 


(6) 


The quantities e* are sampled from a normal distribution with standard deviation S, which we will refer to as the 
disorder strength. This will lead to a unitary ?7' t defining the step of the quantum walk with static disorder. Note 
that time-dependent disorder, i.e. if U' tep depends on the step number, would lead to dephasing of the QW [34]. 
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FIG. 7: Simulation of the IPC implementation of a DTQW in a 2D lattice in non-zero magnetic field with disorder, starting 
at position (X,Y)=(1,1). Probability distribution of a discrete-time QW on a 60x60 disordered lattice with magnetic field after 
50, 100 and 150 steps. The strength of disorder is a = 0.1 and the magnetic held is 0 = 7t/5. The quantum walk propagates 
mainly along the edges of the lattice. The presence of topologically protected edge states that do not localize allows for the 
propagation of the photon to the opposite corner of the lattice. Here we investigate a larger number of steps as compared to 
Figs. 5 and 6, in order to see the propagation to the opposite corner of the lattice, (X,Y) = (60,60). The probability amplitude is 
averaged over 20 disorder realizations. Note that the scales of the plots are different, for better visualization. Total probability 

for each plot: (a) 1, (b) 1, and (c) 0.984. 


Evidence of two-photon edge states 


A quantity that can easily be calculated from the probability distribution of the position of the photons at the 
output of the IPC is the probability that the photons leave the circuit in a waveguide which belongs to the edge of 
the 2D lattice. For the two photon quantum walk, we calculate this probability for a lattice of size 30x30, after 20 
steps, and for variable magnetic flux 0 (see Fig. 8). The two photons are inserted in the circuit in the corner of the 
lattice, position (X,Y) = (1,1), and its nearest neighbour in the x-direction at the position (2,1). The two photons 
are entangled in polarization in a symmetric (antisymmetric) way in order to simulate bosonic (fermionic) statistics 
[28, 34]. We see that, for zero magnetic field, it is very unlikely that the photons leave the circuit by the edge of the 
lattice but with magnetic field this probability increases up to « 15%. It is interesting that the particle statistics does 
not affect much this probability unlike what happens with the average distance between particles shown in the main 
text. 
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FIG. 8: Probability that the two photons are at the edge of the lattice for the quantum walk on the IPC, after 20 steps, 
in a lattice of size 30x30, for different values of the magnetic flux <fi. If the photons’ polarization states are entangled in a 
symmetric (antisymmetric) way, their statistics is bosonic (fermionic). The presence of the magnetic field increases significantly 
the probability that the two photons are at the edge. However, the exchange statistics of the two-photon wavefunction does 

not affect much this quantity. 























